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We explain the dynamics of cold atoms, initially trapped and cooled in a magneto-optic trap, 
in a monochromatic stationary standing electromagnetic wave field. In the large detuning hmit 
the System is modeled as a nonhnear quantum pendulum. We show that wave packet evolution 
of the quantum particle probes parametric regimes in the quantum pendulum which support clas- 
sical period, quantum mechanical revival and super revival phenomena. Interestingly, complete 
' reconstruction in particular parametric regime at quantum revival times is independent of potential 

^ ! height. 
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INTRODUCTION 



Pendulum in quantum mechanics [ij is a subject of great interest when the question comes to explain hindered 
internal rotations in chemistry , quantum features of scattering atoms in quantum optics [1, 0] , perturbation theory 
methods to study weak field effects in quantum mechanics fa-H]' dynamics of Bose-Einstien condensates in optical 
lattices for small nonlinearity[8l| and many other physical systems. Comprehensive study of quantum pendulum has 
, brought to light its various aspects which include the structure of energy spectrum, the time evolution focusing on 
■ the quantum revivals [ol-fl^ and asymptotic Mathieu solutions using algebraic methods [l3| . 
^ \ Super cold atoms in optical standing field is an area of both theoretical and experimental interest and modeled 
O^. as quantum pendulum j3i]- The system is of great importance in classical and quantum domain as an explicit time 
dependence in phase or amplitude makes the classical counterpart chaotic. 

Bose-Einstein condensates in an optical lattice (l^ is a mile stone and an important paradigm to study atomic 
condensate evolution, and its phase transition from Mott insulator to superfluid state |15l42l| . Quantum revivals in the 
system present a profound manifestation of quantum interference. The phenomenon occurs as optical lattice potential 
is not perfectly harmonic and level spacing varies with quantum number. The revivals provide useful Information 
about the coherence time of the atoms in the optical lattice. 

In [ia [TU energy spectrum and time scales encoded in the spectrum are discussed in detail. In this contribution, 
(N • we extend the work foUowing the same approach and explain the evolution of super cold weakly condensed atoms in 
time independent optical lattice. Furthermore, we explain the energy spectrum and eigen states of the system both 
analytically and numerically. The time dynamics of the wave packet also translates the relevance of energy spectrum 
with classical phase space. We probe distribution of energy eigen states of the potential via wavepacket evolution. 
For the simplicity we limit our discussion to deep optical lattice where tunneling does not play significant role. We 
show: (i) An equally spaced local energy spectrum deep in the well which reflects itself in the reconstruction of the 
wavepacket at classical periods; (ii) Beyond this region, the behavior changes as the level spacing modifies itself. The 
nonlinear behavior dominates and controls the dynamics as we go farther from the deep-in-the-well condition; (iii) 
Numerical calculations with enhanced efhciency and accuracy in the presence of analytical relationships explain the 
detailed dynamics of the quantum pendulum and its parametric dependence on anharmonicity. 

In section II, we discuss atomic interaction with optical lattice and obtain quantum pendulum. In section III, 
we work out analytical solutions of the potential. In section IV, simplified quantum pendulum is modeled as series 
expansion of cosine potential and effect of each term of the series is explained which plays the effective role in the 
atomic evolution in quantum pendulum. Results are discussed in section V. 



II. THE MODEL 



We consider super cold two-level atoms interacting with a classical monochromatic standing wave field, E{x, t) — 
ej,[eo cos(A:ia;)e~'"* + c. c.]. Here, w and kL are, respectively, the frequency and the wave number of the laser field, 
and iy is the polarization vector. We assume that the cavity end mirrors, in yz-plane, reflect the incoming light wave 
along a;-axis and, therefore, determine the position of the nodes of the standing light along the axis. In dipole and 
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rotating-wave approximations the atom-field interaction is controUed by the Hamiltonian, 

-ff = ^ + fi^ole) (e| - { d .eyEo cos{kLx) x e'^*a+ + H.c}. (1) 

Here, fiwo is the energy difference between atomic excited state, |e), and ground state, \g), and d indicates atomic 
dipole moment. Moreover, p^; is the center-of-mass momentum of the atom along the cavity axis, M is the atomic 
mass and a± are the raising and lowering operators defined as; (t+ = |e)(f;| and ct_ — \g){e\. Furthermore, we 
have considered condensate densities so low that interparticle interactions are neghgible and single particle approach 
remain apphcable [23j . Thus, we represent the wave function of the system at any time of interaction, t, as, 
\ip {x, t)) = ipg {x, t) \g) + ipe {x, t) e~''^*|e). In the presence of sufficiently large detuning between the atomic transition 
frequency and the field frequency, i. e., Sl = ^^o — i-^, allows us to neglect spontaneous emission. Furthermore, the same 
consideration ahows us to ehminate the excited state adiabaticaUy and leads to effectively describe the evolution of 
the atom within the electromagnetic field in its ground state. Hence the dynamics of an atom, is governed effectively 
by the Hamiltonian, 



H = — -^cosi2kLx). 

2M 8 ^ ^ 

Here, ^eff = ^^/Sl is the effective Rabi frequency, where, = d. eyEo/ h is Rabi frequency. The probability to find 
an atom, therefore, in the excited state is negligible as a consequence of large detuning, the evolution properties of 
the atom in the field are completely determined by the ground state amplitude. It is useful to define the following 
dimensionless quantities t — uit, x — 2kLX, p — {2kL/Muj)px and H — H{Ak\/Mijj'^). We get the dimensionless 
effective Hamiltonian as 

5=Y-^oCOSx, (2) 

where, the parameters Vb = eVl^ff/^'^ is the effective potential depth and e = hk\/2M is the recoil shift. In case 
(5l > O and the laser is tuned red to the atomic transition, we find Vo as positive. Hence, a phase difference between 
the potential and laser intensity shifts the location of the potential minima such that they coincide with the locations 
of the laser intensity maxima. The atom is, therefore, attracted towards the intensity maxima. In the other case, that 
is ^z, < O, the phase difference between the potential and laser intensity disappears and the atom is attracted towards 
intensity minima. The quantized system has another controlling parameter as scaled Planck's constant k = 8e/uj, 
which follows the commutation relation i [p, 5] — k. Here, H effectively defines the dynamics of a quantum particle as 
a quantum pendulum. 

The eigen states for quantum rotor have the spatial periodicity condition (j)n{x) = 0„(a; + 27r), while, eigen states 
for an atom in optical lattice obeys the general Bloch condition (j)n{x) = e*^'^'^0„(a;) where i' indicates the quasi- 
momentum of the eigen state. The general Bloch condition satisfy the spatial boundary condition for quantum rotor, 
when ly has integer values. In contrast, the range of z/ is continuous for optical lattices. In both cases the potential 
only couples the eigen states on a momentum ladder where the ladder spacing is 2hkL. However, for optical lattices, 
there are many momentum ladders independent from each other while for quantum rotor there is only one, centered 

atp=o m. 



III. SCHRODINGER EQUATION FOR QUANTUM PENDULUM AND MATHIEU EQUATION 



The dynamics of the atom interacting with a standing wave field in the large detuning limit is effectively controUed 
by the time independent Schrodinger wave equation, expressed as 

-Y^P- - ^0 = E„^{x), 

where, for simplicity we write ipg — ^ and x ^ x. We rewrite the position variable a; as x = a; — tt, and get Mathieu 
equation, 

^ g^2^'* + K - 2q cos{x)]'iJjn{x) = O, 
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which States that for any given a„ and g, we have a series of solution, ipniz), for the differential equation, labeled by 
index n. Hence, ipni^) defines the eigen function of the system and 

express, respectively, effective potential depth and Mathieu characteristic parameter leading to eigen energies. 

Atom in an optical lattice may observe deep or shallow potential depths corresponding to its energy. Following 
equation |3l we scale effective potential depth by recoil energy of an atom, En = he. The atom in cosine potential 
observes a deep optical potential if effective potential depth Vb is of the order of a few hundred single photon recoil 
energies and temperature of the atom is around recoil temperature. In this case the dynamics of the quantum particle 
in the individual well is independent and one obtains multiple realization of anharmonic oscillators. On the other 
hand when the depth of the potential is just few recoil energies and atom is at about recoil temperatures, it sees a 
shallow potential. In this case, the quantum mechanical effects caused by spatial periodicity of optical lattices such 
as formation of Bloch waves become important. Furthermore the levels are broadened into bands due to resonant 
tunneling between adjacent wells (25| . Tunneling in the low-lying bands is suppressed as the well depth is increased 
and particle motion is dominated by single- well dynamics as it is discussed for large Vg. 

A moderate values of Vg, with an effective Planck's constant of order unity, indicates the deep quantum regime. 
The semiclassical dynamics of the atom in the standing wave field are observed however for large values of Vo , which 
correspond to small values of effective Planck's constant S and here we find several tightly bound energy bands. 
Quantum mechanical effects for small Vb, become important once the atomic de Broglie wavelength ^jA significantly 
exceeds the lattice constant d = i = This givcs the condition 



p2 4^2 ^2 
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For a fixed value of q, there are countably infinite number of solutions, labeled by n. However, only for specific 
characteristic values of the parameter a„, the solutions will be periodic, with periods tt or 2n in the variable x, which 
are denoted by a„, and &„, respectively for the even and odd solutions. Because of the intrinsic parity of the potential, 
the solutions can be characterized as being even, ce„(a;, q) or cosine-like for integral values of n, with n > 0. Whereas 
they are odd, se„(a:, q) or sine-like for integral values of n, with n > 1. Limiting cases i. e.. q = O, q ^ 1 and q ^ 1 
for quantum pendulum are discussed in detail in Reference [Toj . 

Approximate expressions for the characteristic values of a„ and 6„ both in g <C 1 and in g 3> 1 limits are provided 
by references [H, [g^]- For the limit, g <C 1, we find that the a„, 6„ are approximately degenerate for n > 7, that is, 

Un — Un — n -r 2(„2_i) + 32(n2_i)3(„2_4) "t" {'i) 

The above expression is not limited to integral value of n and is a very good approximation when n is of the form, 
771 + i. In case of integral value of n = m, the series holds only up to the terms not involving — m? in the 
denominator. The difference between the characteristic values for even and odd solutions satisfy 

a,i - hn = 0(;;^) as n ^ oo 

In other limiting case, when g S> 1 and the spectrum is oscillator like, we find 

a„«&„+i«-2g + 2sVg-4^-^- , (5) 

where, s = 2ri + f . It has, thus, {n + i)Sa;/i dependence in lower order, which resembles harmonic oscillator energy 
for ujh = 2\A^. Here in the deep optical lattice limit, the band width is defined as [1^ 

n ~ 2^"+^v/^gt + f oxp(-4Vg) , . 

On+1 - , (O) 

which shows that in the deep optical lattice (g ^ I limit) energy bands are realized as degenerate energy levels as 
the band width is negligible. The band structure of optical lattice is shown in figure-1, for large g, near the bottom 
of the lattice, thin band are seen, band width increases and band gap decreases as we move towards the top of the 
lattice potential well. 

As a consequence, we suppress atomic tunneling in deep optical lattice limit. The hoping matrix element. J, explain 
the tunneling between adjacent sites for deep optical lattice [13, [SIHI, viz. 



J = EniVo/ER)-^ exp(-2v%7^). (7) 
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Equations |6l [7] show that the width of the bands corresponds to tunneling of the atom from one lattice site to the 
other. In the limit of deep lattice potentials, this probability will be exponentially small and band width therefore, 
reduces exponentially as a function of lattice potential depth. 

Using algebraic methods, asymptotic even mathieu functions are [ISj given as 



j N _ (Tt"^-6rt-^+llK''-6n)_D„_4 + (4K-4n'')_D„_2-4£)„_|_2-£'n+4 

cen{x,q) — 



^ \^{n^ - 28rJ + 322n^ - 1960n^ + 6769n'' - 1313271^ + 1306871^ - 5040n)D„_8 - («^ - ISn^ + 85n^ 



1024g'8 

-225n^ + 274n^ - 120n)Dn^e + 4(n^ - + 17n^ - \lr? + 6n)D„_4 - (n* + 26n''' - 37^^ + 10n)Ai-2+ 
(-36 - 25n + n^) - 4(n + 2)i?„+4 + + ^^n+s] + A. + 0(4), (8) 

where, in terms of Hermit polynomials -D„(q;) = ^ exp(— ^|-)_ff„(-^)and a is defined as a = 2(7'i cos(a;). The 
normalization factor for ce„(a:, g) up to order O(^) is 

^-~Til+8(^ + 2048^ + ^^^^^^1 + 0(^)J, (9) 

here eigen states are normalized to tt. 



IV. BOUND STATES OF AN OPTICAL LATTICE 



Aroimd the minima of lattice sites harmonic evolution prevails and in the presence of the higher order terms it 
gradually modified to the original potential. Microscopic investigation of the atom-optical system, using term by term 
contribution of cosine potential expansion reveals the dominant role of the system's particular parametric regime in 
the formation of eigen-states and eigen energies. This leads to simplified analytical solutions around the potential 
minima in the system, as discussed below. 



A. Harmonic Oscillator Like Limit. 



In the deep lattice limit, the potential near the minima can be approximated as quadratic. Thus the particle placed 
near the minima of the cosine potential, experiences a harmonic potential. The energy in this regime is obtained as, 
= (2n + l)Sv^ - Vb, which can be identified in equation[5l by ignoring square and higher powers in s. The 
eigen states of quadratic potential are given as 0„(a;) — y^ ^„^^^^ iJ„(/3a:) exp( ), where, 7?„(/3x) are Hermite 

polynomials and /3 = (^^^)^ = \/2'?^- 

The time evolution of the particle, initially in state "(/"(s;, 0), is obtained by time evolution operator J7, such that, 

oo 

^{x,t) — U'ip{x,0) — ^ c„0„(a;) exp(— where, £"„ and 4>n{x) are energy eigen values and eigen states corre- 

Tl = 

sponding to quantum number, n. The probability amplitudes Cn are defined as, {4)n(x)\ip(x, 0)). The quantum particle 
wave packet in optical potential, narrowly peaked around a mean quantum number n, displays quantum recurrences 
at different time scales, defined as T(j) = f , where, E^^ denotes the j*^ derivative of En with respect to n. 

The time scale, 7(1), is termed as classical period as it provides a time at which the particle completes its evolution 
following the classical trajectory and reshapes itself. Whereas, at T(2) the particle reshapes itself as a consequence 
of quantum interference in a nonlinear energy spectrum, which is purely a quantum phenomenon and thus named as 
quantum revival time. In the parametric regime of a changing nonlinearity with respect to quantum number, n, we 
find, the super revival time T(3) for the quantum particle [29[. 

We study the time evolution of the material wave packet using square of the autocorrelation function [30l |. 

I^WP = E + 2 E ICnHcmP COs[(i?„ - E„,)-]. 

n— O n^m 
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In the present parametric regime the square of the autocorrelation function is written as 

oo oo 

I^Wl' = I]|cn|^ + 2^ \cr,Wm\^cos{{n-m)2^ot), (10) 

n— o n^ra 

oo 

where, ^ |cn|^ is independent of time, and defines the interference free, averaged value of The most dominant 

n=0 

contribution to the |A(i)p comes from n — ra — \ the second part at the right side of eauation llOl Other terms with, 
m — n > 1, have negUgible role because their osciUation frequency is an integral niuhiple of fundamental frequency 
and are averaged out to zero. For the reason the square of the autocorrelation function in the present case 
oscillates following a cosine law with a frequency 2^/Vq which leads to a classical time period T^^[''' = as shown in 

figure 2. Here, zero in the superscript of Tj.^^ defines the system's classical period in the absence of perturbation. At 
the integral multiples of classical period, is unity, whereas, at time, which is odd integral multiple of the half 

of the classical period 

oo oo 

\mf = E + 2 E Icni'lcml' cos[(n - m)n]. 

n=0 n^m 

Here cos[(n — m)TT] has altcrnatively values +1 and —1, whcn n — m is an even or odd respectively. Thus after the 
cancellation of positive terms with the negative ones, second summation reduces to minimum value and |^(i)p attains 
its minima. Here, in the eigen states expansion 

tp{x,t) = exp[~i{2^/% - -^)t]^Cn(t>n{x)exp{-i2n^/\^t), (11) 

n=0 

it is notable that the eigen states 0„(a;) have parity (—I)"- In case of even parity, only the even terms C2n are 
nonvanishing and n-dependent exponent factor oscillates two times faster than the general case. At half of classical 
period the wave packet reappears towards other turning point of the well. In case the wave packet is initially placed 
at the center of the cosine well, it reappears at half classical period at the same position but in opposite direction. In 
this case we see classical revivals of the initial atomic wavepacket and quantum revivals take place at infinitely long 
time. Hence, we find revival of the atomic wavepacket after each classical period only. Experimentally we may realize 
the situation by placing very few recoil energy atoms deep in the cosine potential well. This reveals the Information 
about the level spacing around the bottom of the cosine potential. Interestingly we find an equal spacing between the 
energy levels from figure 1, for large q and small value of n. The spatiotemporal behavior of the wave packet in the 
quadratic potential, as shown in figure 2, confirms the above results. 



B. Quartic Oscillator Limit. 

Beyond harmonic oscillator limit, we find oscillator with nonlinearity and energy level spacing different from a 
constant value. The correction to the energy of harmonic oscillator comes from the first order perturbation (see 
appendix for energy corrections), that is, E'i^^' = — ^(2n^ + 2n + 1), which again can be identified in equation [5] by 
ignoring cubic and higher order powers in n. The atoms with little higher energy, which is equivalent to several recoil 
energies observe another time scale in which it reconstructs itself beyond classical period, i. e. quantum revival time. 
The behavior of auto-correlation function for the wavepacket exactly placed in this region, where only first order 
correction is sufficicnt, is shown in figure 3. We see that the wave packet displays revivals at quantum revival time. 
Thus, little above from the bottom of the of an optical lattice, the wave packet sees variations in energy level spacing in 
a non- linear fashion. The classical time is modified as T^^P = a'^^'T^^^-*, where a'-^' = 1 + and classical periodicity 
for a particle in present situation is related to potential height and mean quantum number, here, s = 2n + 1. As n 
increases, classical revival time also increases, and the ratio, is always less than unity in the region where first 

correction is sufficient. The quantum revival time, Trll = is independent of mean quantum number n, whereas 
super revival time in this case is infinite. 

The eigen states of quartic oscillator due to first order perturbation are given as 0^(a;) — 4'n{x) + (x), where, 
ix) is first order correction to the harmonic oscillator wave function and is defined as 



(12) 
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Here, Di, rji, 772, Vs ^-^d 774 are constants and defined in appendix. 

In figure 4 eigen states of quadratic, quartic, sixtic and octic oscillators are mapped on numerically obtained eigen 
States of cosine potential. From this figure we see that for Vb — 10 and S = 0.5, first order correction to the eigen 
states of unperturbed system matciies with harmonic oscillator eigen states up to 71 = 3, and mapping of quartic 
oscillator with exact solution is much improved compare to harmonic oscillator. Siniilarly mapping of sixtic oscillator 
is better than quartic oscillator and is quite good for octic oscillator for all bound band. In this case eight bands exist 
in side the potential. From figure 1 we note that by increasing q number of bound bands can be increased. 

The square of autocorrelation funetion in this case takes the form 

00 00 ^ 

\A{t)\^ = ^|c„|4 + 2 ^ |c„nc™pcos[(n-m)2v/T^t+(m-n)(n + m+l)-t]. (13) 

n— o n^m 

Here, nonlinear dependence of energy eigen values on quantum number, n, makes the argument of cosine funetion non- 
linear, as appears in the last term of the above expression. The nonlinear term (m— ri)(n+m+l) removes a degeneracy 
present for the harmonic case between the cosine waves corresponding to nearest neighbouring off diagonal terms and 
beyond. Hence, there overall evolution display a gradual decoherence leading to coUapse which latter transforms in 
revival as the decoherence in waves disappears. The values of in this regime at Trev — ^ simplifies as 

00 00 

I^Wl' = I]|cn|'' + 2^ |c„|2|c™|2cos[16^y^(n-m)]. (14) 

n=0 n^m 

The equation[T3]displays that |^(t)P = 1, when ^ is an integral multiple of i and also |^(t)P is unity at half of the 
revival time if ^Jq is an integral multiple oi\. In case the ^Jq is not an integral multiple of | , the wave packet revival 
occurs little earlier than the revival time Trev — ^ and also |^(i)P approaches to unity little earlier than Here 
reproduction of wave packet at is out of phase by tt, i. e. at that time all the waves are moving exactly in opposite 
directions as initially they were. But at T^-^v they are all moving in same direction and each wave is in phase not only 
with there initial states but also with each other. 

Figure 5 shows spatiotemporal behavior of the wave packet in cosine potential. We see that wave packet spreads 
and oscillates in the cosine well, and after some time the original wavepacket is divided into sub-wave-packets figure 
5 (b). At a quantum revival time these sub-wave-packets constructively interfere and wave packet gains the original 
shape at the same initial position figure 5 (c). At quantum revival time, same classical pattern is seen as at the start 
of the time evolution shown in figure 5 (a). 

C. Sixtic Oscillator Limit: Existence of Super Revivals and Beyond 

Higher order nonlinearity in the energy spectrum of the quantum pendulum show up beyond quartic limit. In the 
presence of the second correction, the energy is modified by the term E^n^ = — (2k +3n^+3w+i) ^ gg^^j-^^j correction 

to the energy modifies the time scales. The classical time period is now = a^'^^T^^ ^ whereas the quantum revival 
time is modified as Trew = \^'^^^\TrVv- Here, a'-^' = a^^^ + -^^^r^, ^'^^^ — ^ — 1 are constants. In addition system 

shows another time scale, i. e. super revival time Tjpr — ^^^3-^, at which reconstruction of original wave packet takes 
place. The other quantum revival times occur at infinity. 

The energy eigen states in this regime are given as ^i"'' — 0« + + ^i^'''' -f 4>n'°'\ where, is perturbation 

in eigen states due to the H^^"^ term and ^i^'^'' is second order perturbation caused by i?^*) term [T3|. The expressions 

,(l,b) J ,(2, a) 

(pn and (pn are given as 
and 

The constants D2, Dg, S j' s and Xj's are calculated in appendix, where, j takes integral values. 

Again in this region classical time increases by increasing n and increases faster than as it was in quartic limit. 
However, in the present regime, the quantum revival time is not constant and decreases as n increases. The n 
dependence of quantum revival time is shown in figure 6. The super revival time is independent of mean quantum 
number. It is directly proportional to square root of potential height and inversely proportional to the square of 
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scaled Planck's constant. The temporal behavior of an atom in optical potential placed in this rcgimc shows three 
time scales: classical periods enveloped in quantum revivals and quantum revivals envclopcd in super revivals as 
shown in figure 7. After each super revival time the atomic wave packet evolution repcats itsclf. 

Similarly third correction in energy, modify the cncrgy by the factor En"^ ~ ~% +Wn +^16« +iiw+3) rp-j^^ time 
scales in this case are T^f = a^^^rj,"^ and T^^^ = IP^^^T^el Where a(3) = a(2) + ^^2!±^ and /3(2) = /3(1) + 

and super revival time is Tspl = {jlTepl where 7 = — 1. Furthermore, the super quartic revival time, T4, is 
independent of n. The higher order corrections in energy show that other time scales do exist in the system, but their 
times of recurrence are too large to consider them finite. 

In the prcscncc of third correction to cncrgy, classical time incrcascs as h incrcascs but increases littlc fastcr than in 
the cases of quartic and sixtic corrections, whereas, the quantum revival time decreases faster as n increases compared 
to the case of sixtic correction. The super revival time is not constant but decreases as n increases. 

Energy corrections incrcasc anharmonicity in the system. Wc discussed that the first order correction to harmonic 
potential energy led to the quantum revivals, the second order energy correction led to super revivals and fourth 
correction led to quartic revival time. Comparison of revival times for different energy corrections is shown in figure 
6. For first order energy correction quantum revival time is constant, but for higher order corrections, revival time 
may decrease with increasing mean quantum number of the wavepacket. 

In figure 4 we show the projection of numerically calculated eigen states of cosine potential on the eigen states of 
quadratic, quartic (first correction to quadratic potential), sixtic (second correction to quadratic potential) and octic 
(third correction to quadratic potential) potentials. We note that the eigen states of all above mentioned potentials 
match to the cigcn states of quadratic potential ncar the bottom of lattice potential. For littlc large quantum numbcrs, 
we see that projection of quadratic potential falls sharply and improves with higher order potentials. This correction 
is quite good for octic potential when q = 40 justifying q » 1 condition. Eigen states and eigen energies of quartic, 
sixtic and octic potentials are analytically calculated using the method given in appendix. 

V. DISCUSSION 

In this paper we have extended the understanding of eigen energy levels and cigcn states in deep optical lattice, both 
analytically and numerically. We note that solutions obtained through perturbation theory and Mathieu solution are 
showing similar results and are in very good agreement with exact numerical solutions. The energy levels are equally 
spaced near the bottom and by increasing the lattice potential depth, the number of cqually spaced energy levels can 
be increased. A wave packet placed in this region revives after each classical period. From figure 4, it is also noted 
that all potentials discussed have equally spaced eigen levels at the bottom of the potential well as their mapping 
with harmonic oscillator is unity. Beyond linear rcgimc, energy dcpcndence is quadratic and any wavepacket evolved 
in this region shows complete quantum revivals enveloped by classical revivals. Interestingly in this regime quantum 
revival time {Trll = ^) is independent of potential height, however has inverse proportionality with effective Planck's 
constant, S. We show that for dcep optical lattice, there is a region in which revivals are independent of lattice depth 
and this region expands with the increase in lattice depth, however beyond this region quantum revival time is no 
longer constant but decreases as n increases. The higher order time scale, super revival time Tgpl = ^^^^^ exists in 
this region and is independent of n. Again this region cxpands as potential height is increased but super revival time 
in this region is directly proportional to square root of potential height. 

Above this region other time scales also exist where super revival time is n dependent and decreases as n increases 
but these time scales are too long to consider them finite. 
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VII. APPENDIX: SOLUTION OF AN ARBITRARY POTENTIAL 

An arbitrary potential U{r), around its minima can be solved by taking its Taylor's expansion ^31j . that is 

U {r) = t/(r„) + G(i)(r-r„) + G(2)(r-r™)2 + G(3)(r-r™)3 + ...., (15) 

where, G(^) = {i^)-^dW{ T = rm)/dr^ , and j is an integer. The value of G^-'^ for odd j is zero as potential is 
cos(x) and it is calculated at the potential minima r = r^. Thus in the presence of weak nonlinearity the term 
G(6) «G(4) «G(2). 

In our analysis we consider the expansion up to second order term as unperturbed Hamiltonian, Hq. The eigen 
functions and eigen energies of this Hamiltonian are those of harmonic oscillator. The efFect of the higher order terms 
in Taylor's expansion is discussed as perturbation to eigen energies and eigen functions of the harmonic oscillator. 
We express the effective Hamiltonian governing the dynamics of atom around the potential minima as 

^o = y + f/(r™) + G(2)(r-r„0^ (16) 



The eigen functions and eigen energies of the harmonic potential are 4>n{x) = ^ ^„^^^ -ff„(/3a;) exp( ), and 

En'^ = S-yVo^(2n + 1) + U{rm): where, Hn{(3x) are Hermite polynomials. 

The first order correction and second order correction to energy is quite well known. The first order correction 
to the energy of quadratic potential is E'i^'"'' = (0„|i/^P-' |0„), and the second order correction is obtained from the 
relation E^'^^ = X]m=^n ^^''^"m' ^^t»?^^ • "^^^ third order correction is determined by 

— -£'771 

_ ^^ (n|g(rt|j-)(^-|g(p)|^)(;|H(P)|n) ^ (^n\mP)\j){j\H'^P)\n) 



i^nl^n {^n - Ej )(En - ii; j j-^„ [(En ) - [Ej ) ) 

where, p = 3, 4, 5, 6, • 

The presence of the perturbation term H^^\ leads to the Hamiltonian 

H = Ha + H'^^\ (17) 
Here, the first order correction to eigen functions, due to the ff term, is obtained as 

= L)i(?7i0„_4 + 7/20n-2 - m4>n+2 " r]4.(t)n+4) , (18) 

where, Di = gW(-J^)2 i 



m = \/nin-l){n-2){n^3)/A, 
772 = (2n-l)v/n(n-l), 



7^3 = (2n + 3)v/(n + l)(n + 2), 



and 774 = ^(n + l)(n + 2)(n + 3){n + 4)/4. 
Now the eigen functions in the presence of first order perturbation, due to the correction H^'^\ is given as 

and second order correction due to iJ*^*^ term is 

4>n'''^ = -D2[^10„-8 + S2<f'n-6 + '53(/)„_4 + (540ri-2 + '55(?!)„+2 + '56(/>„+4 + (570n+6 + ^80n+8- (19) 
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y/n (n - 1) (n - 2) (n - 3) (n - 4) (n - 5) (n - 6) (n - 7) 



^1 



32 



^2 = - 1) (n - 2) (n - 3) (n - 4) (n - 5), 

^3 = (2n^ - 9n + 7) v/n(nrny(rr^2Hrr^, 

(56n3 - 228n2 + 214n - 146) ^- 

^4 = ^ g '-^Jn(n-l), 

(56n3 + 396n2 + 838n + 645) , — 

^5 = ^ ^ '-^{n+l){n + 2), 

(31n2 + 197n + 258) r-. —, — —, 

^6 = ^ T- V(«+l) (" + 2) (n + 3) (n + 4), 

lo 

5j = 2^ ^ + 1) (n + 2) (n + 3) (n + 4) (n + 5) (n + 6), 
and ^8 = '^"+^) ^"+^^ ^"+^) ^"+^^ ^"+^) ^"+'^^ ^"+^ 

Hence, foUowing the same procedure, the first order correction due to i?^^^ term changes the Hamiltonian of the 

system as 

H = Ho + H^"^^ + H^'^l (20) 
Hence the corrected eigen function in presence of correction due to H^^^ and H^^^ terms appear as 

where, 

Here, D, = G^^) {^f 



Xi = 6Vn(n - l)(n - 2)(n - 3)(n - 4)(n - 5), 

X2 



^(2n - 3)V«(n-l)(n-2)(n-3), 



X3 = yln^ - n+ 1)V"(" - 1), 

X4 = Y(n2+3n + 3)V(n+l)(n + 2), 

X, = \{2n + 5) V(n + 1) (n + 2) (n + 3) (n + 4), 

X6 = 2V(n + 1) (n + 2) (n + 3) (n + 4) (n + 5) (n + 6). 

Close to the minima of potential, we can find the eigen energies up to a considerable accuracy by using perturbation 
thcory. The Icading correction comes from i?^*) using first order and second order perturbation theory respectively. 
The result is as under: 

eI^^ = (a2«^ + ain + ao)'kLJh, 

where, «o = 3£'a, = 6Da, = 6D„, Cb=iG^^''f{j^f^ and = G^^Hj^f)^- At next order, the first order 
perturbation of il^^^and second order perturbation of H^'^^ contribute. The result can be written as: 

where, po = 37« - 21Ji, = 87« - 59 J^, p2 = 67„ - SIJ^, P3 = - 34 J^, and 7„ = 5G(6)(_i=)3_i_, = 
2(G(^^)^(i^)^(^)^- 
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At the ncxt highcr order, wc need to evaluate three contributions; H^^^ in first order, H^^^ and H^^^ in second order 
and H^^^in third order. The energy expression is 

En'' = (74«'' + lan^ + 72«^ + 71" + 7o)Swft,, 

where, 70 = 3X - 12Y - lUZ, 71 = 8X - 35Y - 347Z, 72 = 10X - 46Y - 472Z, 73 = 4X - 22Y - 250Z, 
74 = 2X - IIF - 125Z, and X = 35G(8)(^)4^, y = 30G(6)gW z = 48(GW)3(^)6(_i_)2. 

Now the energy of the system is 

= E^^) + E^:) + E^:) + E(:\ or 
En = {KiU^ + K^n^ + K2'n?' + K.\n + Ko)TiuJh + U{rm), (21) 

where, «o = «o + /3o + 7o + 5, k^i = + /3i + 71 + 1, «2 = 02 + /32 + 72, «3 = /^s + 73, and K4 = 74- 
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FIG. 1: Characteristic values for the Mathieu equation, a2m (for even solutions, blue (grey) curves) and b2m (for odd solutions, 
red (dark) curves, mostly overridden by blue curves) versus q for the quantum pendulum: The dotted lines correspond to 
characteristic parameter = ±2g i. e. E = ±Vb- Note that for a given value of q where q ^ a-n, the gap between the lowest 
energy state (lowest solid curve) is roughly one-half of the spacing between solid and dashed curves, corresponding to the 
zero-point energy in the oscillator limit. 
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FIG. 2: Time evolution of particle wave packet placed at the bottom of the cosine potential. The dimensions of the wave packet 
are 'k = 0.5, Ap = 0.5 with Vb = 10- We show autocorrelation function vs time (right side) and spatiotemporal behavior of 
the material wavepacket (left side). The wave packet see equally spaced energy levels and rebuilds after every classical period. 
Analytically calculated value of classical period and numerical results show an excellent agreement. 
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FIG. 3: Autocorrelation function for a particle undergoing quantum revival evolution in time. The parameters are same as in 
flgure 2. The wave packet was placed close to the bottom of the potential well in the regime where first order correction is 
sufflcient, it observes quantum revivals after many classical periods. 
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FIG. 4: comparison of cosine potential with the simplified potentials is made by calculating the projection, S, of the eigen states 
of the cosine potential on the eigen states of the simphfied potentials . For first few quantum numbers, the cosine potential very 
much resembles to the harmonic potential, however, for higher quantum numbers the higher order corrections to the harmonic 
oscillator are needed to make the resemblance. 




FIG. 5: The time evolution of a localized wave packet in a cosine potential is displayed for the same parameters as in figure 
2. For short times the initial wavepacket shows classical revivals (a), but starts to display sub-wave-packets in its long time 
dynamics (b), which constructively interfere at quantum revival time Trev (c). Our analytical and numerical results are in very 
good agreement. 
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FIG. 6: The quantum revival time vs mean quantum number is shown for simplified potentials. In the presence of only first 
order correction, i. e. for quartic potential, the quantum revival time is a constant. For higher order corrections it decreases 
with increasing mean quantum number. 
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FIG. 7: The wave packet dynamics in sixtic potential displays three time scales, the classical periods (making the dense region), 
the quantum revival times (making the peaks in the dense region) and the super revival times (making the peaks in the envelop 
of the quantum revivals). The parametric values are the same as in figure 2. 



